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Abstract—This paper first provides an intuitive description 
of the spectral representations for the transmittance between 
two apertures. Then, two spectral representations of the trans- 
mittance between OAM antennas are compared. The first one 
makes use of the Modified Zernike function and has already been 
described in a previous publication. The second one is based on 
the Fourier-Bessel decomposition of the current distribution on 
the apertures. Both families of expansion functions are orthogonal 
on a disk and admit a closed form spectrum. In both cases, the 
transmittance in the far-field can be approximated in closed form 
if the current distribution is expressed in the basis. However, it 
is observed that the far-field behavior is fundamentally different, 
highlighting the complementarity between the two bases. 


Index Terms—OAM, aperture antennas, transmittance, multi- 
plexing, spectral domain formulation. 


I. INTRODUCTION 


Rthogonal Angular Momentum (OAM) communication 

[1] has been presented a few years ago as a possible new 
way to multiplex signals within a given frequency band. In a 
nutshell, each channel corresponds to a specific phase roll by a 
multiple @ of 27 versus azimuth. If transmitting and receiving 
antennas can support several modes, independent signals can 
be transmitted within the same frequency band, as is already 
exploited for several decades for optical fiber communication 
[2]. A first demonstration in free-space RF communications 
has been described in [1]. It is interesting to notice that the 
generation and, above all, reception of such modes requires an 
aperture with substantial size. Moreover, higher order modes 
decay faster with distance, which limits the process to the 
near-to-intermediate field [3], [4]. 

The space-spectrum domain provides a relatively synthetic 
way to look at this problem. In [5], a general formula has been 
provided for the transmittance H between two apertures in 
the form of a mutual impedance (see Fig. 1). It is generalized 
below to directly include the vectorial aspect of fields, through 
the use of the dyadic Green’s function G(r, r’) for free-space 


[6, eq. (7.2)]: 
— 4k > 
l-l S,. 
where K, and K, are respectively the equivalent vector current 
distributions on the transmitting and receiving antennas, as- 
suming that both antennas are active. /; is the current flowing 


at the transmitting antenna port, 77 is the free-space impedance 
and k is the free-space wavenumber. The link between the 
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transmittance, and the power transfer has been proven in 
Appendix C. 

Equation (1) can be interpreted by discretizing the integrals 
over transmitting and receiving apertures, such that the current 
distributions can be replaced with virtual antenna arrays, 
characterized by amplitude distributions K, and K r- If the 
arrays are placed very far from each other, at the level of the 
receiving array, the fields from the transmitting array strongly 
resemble a plane wave, which is the only communication 
channel between the two apertures. When the curvature of 
the phase front at the level of the receiving antenna be- 
comes non-negligible, this single-channel model is no longer 
valid. Indeed, the curved phase fronts generated by different 
transmitting antennas appear as significantly different at the 
receiving array level (dashed in Fig. 1). Hence, after proper 
beamforming, represented by distribution K,., signals from 
different transmitting antennas can be distinguished. This is 
also true when individual transmitting antennas are replaced 
by the whole array, with different distributions K +. Proper 
sets of distributions K, and K,. create several independent 
channels (see Fig. 1) and this function is also exploited in 
free-space MIMO communications [7], [8], [4]. The number of 
independent channels increases as the two arrays (or apertures) 
are approached. 





Fig. 1. Illustration of the different channels created between two apertures 


It is convenient to provide a spectral-domain description of 
the field distributions, i.e. along spatial spectral coordinates 
(kx, ky). Coming back to continuous field distributions, as the 
apertures are approached, a wider set of samples in (kg, ky) 
plane can be exploited to create independent channels. Each 
sample corresponds to a different harmonic in space domain. A 
change of basis can now be operated: instead of decomposing 
the fields in harmonics versus x and y coordinates, one can 
decompose them into a family of functions forming a complete 
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set over the disk. If, besides, those functions are orthogonal 
to each other, the coefficients multiplying each element of 
the basis can be simply obtained by projection. Two such 
families of functions exist. The first one corresponds to the 
Zernike functions [9], [10] and has already been exploited in 
[5] for the study of the transmittance between OAM antennas. 
The second family of such functions is the Fourier-Bessel 
family [11]. As will be seen later, although the derivation is 
similar, the convergence of the transmittance for both families 
of functions — versus number of modes and versus distance 
— have fundamentally different behaviors. Depending on the 
application one may thus prefer one or the other. 

The remainder of this paper consists of deriving analytical 
expressions for the transmittance between antennas whose 
distributions are described by Fourier-Bessel functions. It 
will be assumed that there is no multiple-scattering taking 
place between the apertures, i.e. the fields in and around 
the transmitting antenna are not significantly affected by the 
presence of the receiving antenna. Asymptotic expressions will 
be provided for large distances and large function indices. The 
paper will be concluded by comparing the convergence of the 
transmittance versus the index of either Zernike or Fourier- 
Bessel series and by highlighting an essential difference in 
the intermediate-to-large distances range. 


Il. MATHEMATICAL FORMULATION 


To calculate the transmittance (1), it is more convenient to 
project the currents onto two orthogonal axes. In this paper, 
we use a Cartesian system of coordinates, with unit vectors 
(âx, ây) (see Fig. 2). The transmittance can then be described 
as a matrix standing for the different polarization pairs. As 
proven in [5], the cross polar transmittance is equal to O. 
We will therefore limit ourselves to the co-polar components 
(meaning the diagonal elements of the transmittance matrix). 


A. Fourier-Bessel functions 


The Bessel functions form a class of orthogonal basis func- 
tions (see Appendix A), defined on a unit disk (0 < po < 1) 
as follows [11]: 


Fma po) = JalAnmpo) (2) 


where Anm is the m-th positive zero of J,(x). This family of 
Bessel functions admits a closed-form spectrum (see Appendix 
A). 

The azimuthal behavior is taken into account by adding 
the factor e~J"%, where a is the azimuthal coordinate. The 
Fourier-Bessel functions are then written as: 


Rn (2a) =Faa(®) MO 


where a is the disk radius. 

The resulting family of Fourier-Bessel functions is rigorously 
orthogonal on a disk and admits a closed form Fourier trans- 
form: 


Ra n( 0s) = ~anjrerins | omdn-s om) (9) 


an _ OF 
(4) 


where (@ and ¢ are the spectral variables in cylindrical coor- 

dinates (ky = 8 cos ġ and ky = £ sin @). 

The scalar current distributions (obtained after projection in 

the orthogonal basis) at the transmitting and receiving antennas 

can then be expanded into the Fourier-Bessel basis as follows: 
N 


M / 
K (r”) = `“ > Ve Elma (Za) (5) 


n=—N m=l1 


N M 
k= S Erea a) O 


n=—N m=1 


The systems of coordinates (p’, a’) and (p, a) are represented 
in Fig. 2. a and b are respectively the radius of the transmitting 








Fig. 2. Geometrical configuration: two apertures along with the related 
geometrical parameters. 


and receiving antennas. The expansion coefficients 7, and 
Yr,,, can be calculated by a simple projection: 


Vtnm = a f K (r) RY f a’) p'dp’ da 
(7) 


where K is a normalization factor given by: 
K(m, n) =- datil Anm) Jari Anm) a" (8) 


A similar expression is obtained for yr„„. As an illustrative 
example, Fig. 3(a) and Fig. 3(b) show respectively the absolute 
value of the current distribution and that of the x-directed 
current distribution on a circular modulated metasurface an- 
tenna radiating a broadside pencil beam as analyzed in [17], 
[18]. The real parts of the modes (m,n) = (1,0) and 
(m,n) = (2,1) are represented in Fig. 4(a) and Fig. 4(b) 
respectively. 

Finally, the co-polar transmittance between the two OAM 
antennas, assuming the same azimuthal behavior (n index) 
of the current distribution can be decomposed in all pairs of 
transmittances as follows: 


H= Y ie an anin (9) 


where each contribution of the co-polar transmittance is given 
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Fig. 3. a) Absolute value of the surface current distribution on an anisotropic 
MTS analyzed in [17], [18] and used to generate a right-handed circularly 
polarized (RHCP) broadside pencil beam. b) Absolute value of its x-directed 
component. 
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Fig. 4. Real parts of two modes of the antenna shown in Fig. 3: a) (m=1, 
n=0) and b) (m=2, n=1), where (m,n) has been defined in expression (3). 


by: 


= Hanne =C fS [f Bnn (La) GEP) Po 
Sr St b 
Rmn (2. a’ 
a 


In this equation, C = —j4kn/(I,I;) is the constant outside 
the integral in (1) and G(r, 7") is the scalar free-space Green’s 
function. Ps is a correction factor that stems from Green’s 
dyadic and takes into account the polarization. It is interesting 
to mention that expression (10) resembles the well known 
reaction integral appearing in the classical Method of Moments 
formulation [17]. 


B. Plane waves expansion of the Green’s Function 


The scalar free-space Green’s function can be written as a 
spectrum of plane waves [12], [13], [14]: 





G sah eg e-iklir-r']] 
OP) aT 
1 eI Bp cos(a’—¢) »—JjBp cos(a—$) »—jkzd 
— 2 JJ — lo PD 
(27) 2j k, 
(11) 
where k = 2r/A is the free-space wave number, 8 = 


k2 +k? and ką = yk? — 8?. To take into account the 
polarization effect (treated through Ps in (10)), we need 
to eliminate the contribution of the source parallel to k, 
before projecting the field on the polarized current at the 
receiving antenna. Indeed, (11) corresponds to a plane wave 
decomposition and for each plane wave, the electric field must 
be orthogonal to the wave vector. It has been proven in [5] 
that this is equivalent to multiplying each plane wave of the 
expansion of the scalar free-space Green’s function by the 


factor: 
o (&.RE\ | 


where €; and €, are respectively the transmitted and received 
current polarizations. After inserting (12) into (11) and (11) 
into (10), and using the following expression: 





(12) 


27 
J oina iB cos(a’ —$) Ja adr In(Bp') eI j” (13) 
0 


we obtain (after integrating along @) for the co-polar transmit- 
tance: 


/ —jkz 
n=c= | | Ja (Bp) Emn (2) o do £ 
J B Jo! a k 


J TABOE (2) p dp B (1 = =) dB 
Í (14) 





Using the closed form Hankel transform of the Bessel func- 
tions (see Appendix A), the transmittance can be rewritten as: 


n _ Crab” [ 
J B 


e-ikzd 


Fw (Ba) Fa (Bb) =— 6 
r=0 z 


BP 
(1-8) 


To avoid the pole at 6 = k and to improve the integration 
convergence, we choose as in [5] a very simple linear contour 





(15) 
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given by 3 = (6, (1+ 7), where 6, is the real part of 6. The no effect on the slope with the Bessel decomposition. It can 


final expression of the transmittance is the following: be proven that the dependence of H versus m in the far-field 
CSS [9 ~ 7 B2 is proportional to 1/m° for high values of m (see Appendix 
fe J F? (Ba) F? (8b) (1 — if] B). This is a fundamental difference with the Zernike family, 
= Pr-O l i (16) Where n and m both affect the slope of H versus distance in 
e d a “log-log” scale (see Fig. 8). 
Pr Pi 


II. NUMERICAL RESULTS 


The chosen parameters for the examples below are: D = 
2a = 2b = 40 cm, the wavelength \ = 1 cm and the distance 
d is normalized w.r.t. the far-field distance d fa, = 2D?/). 
Fig. 5 shows the evolution of the transmittance for different í ~ 
values of the n index in a “log-log” scale. We can observe the fett Trien 
same trend as with the Zernike basis (see Fig. 6). The n index 
controls the slope of the transmittance in the far-field. 


\ = = “fey, 





d/d 


far 


Fig. 7. Transmittance versus distance w.r.t. far-field distance. n = 0 and m’= 
m = 1...4 (Fourier-Bessel basis). 
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Fig. 5. Transmittance versus distance w.r.t. far-field distance. n = 0...4 and 
m’ = m = 1 (Pourier-Bessel basis). 
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Fig. 8. Transmittance versus distance w.r.t. far-field distance. n = 0 and m’= 
m = 1...4 (Zernike basis). 


IV. CLOSED FORM EXPRESSION OF THE TRANSMITTANCE 


Under the far-field approximation, the transmittance expres- 
i ó i sion can be well approximated with a closed form expression. 
10 10 10 First, let us replace in (16) the spectrum of the Bessel 
d/d; functions (22). 
2C StSp f Aniw n= Ani A (Ga) 
p 


Fig. 6. Transmittance versus distance w.r.t. far-field distance. n = 0...4 and H = —— X 2 — ( Ba)? 
r=0 nm’ 


m’ = m = Q0 (Zernike basis). a 
Nina) Ja 00) e Ikea ae 
However, as can be seen in Fig. 7, changing the value of m has (Anm)2 — (Bb)? k. Pr (1-3 P dpr 





(17) 
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For a large distance d, the e~/*% factor can be approximated 
as e~Jk4 @—Brd/ which explains why the integrand converges 
very fast. As a consequence, the terms Ga and 8b appearing 
in the denominator can be neglected, in comparison with 
Anm’ and Anm respectively and we can assume that the k, 
factor appearing in the denominator is almost equal to k. 
Expression (17), can then be approximated as 


s CA enh idt f 
p 


H In (Ba) In(Bb) 


TR Naa Anm =0 


—jkd ~—B2d/k Be 
€ e Br 1 “Jp dpr 
(18) 


Next, for the same reason as before, the Bessel function can 
be decomposed using its first order approximation (J„n(£) ~ 
(x/2)"/n!; n > 0). Finally, using the following identity: 


0° ; Le 
J oe Brd/k gl+2ngg — > (5) (19) 
Br=0 


the transmittance takes the following asymptotic form: 


festenan Aee] gana 


nm’ Nea n! And 


(Sr) (wena) 


is the correction factor taking into 





H & ACS; S, | 
(20) 





Expression € — tia} 
account the polarization effect stemming from (12). We can 
observe that this factor tends to 1 when the two antennas are 
far from each other, meaning that in this case, we can assume 
that the polarization of the transmitted field is homogeneous 
over the reception aperture. The expression between brackets 
shows the m-dependency of the transmittance in the far-field. 
It is proven in Appendix B that this dependency is proportional 
to 1/m? for n < m and m =m’ !. Fig. 9 shows the trend of 
the transmittance, assuming the obtained far-field closed-form 
approximation. 

We can observe a very good agreement down to about half 
the far-field limit. A similar approximation has been obtained 
in [5] for the Zernike family (see Fig. 10). In the latter case, the 
propagation exponent was —(1+n+2m) instead of —(1 +n). 


V. CONCLUSION 


Two formulations have been proposed to study the trans- 
mittance between OAM antennas. In both cases, for large 
distances, an asymptotic expression has been provided. Both 
formulations cater for the radial variation of fields across 
the apertures. This has been represented using two different 
families of functions orthogonal over the disk, i.e. Zernike 
and Fourier-Bessel functions. Roughly speaking, the number 
of zero-crossings of the functions along the radius is given 
by the m index in both cases. However, it is interesting 
to notice that the dependence of the transmittance on m is 
completely different. When Zernike functions are used, index 
m appears (with a factor 2) in the exponent of the decay 
rate. When Fourier-Bessel functions are used, m does not 


lFor m ~ m’, we have m~3/? and (m’)~3/? dependences. 
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Fig. 9. Transmittance versus distance normalized w.r.t. far-field limit. m’ = 
2,m = 1 and n = 0...4 (top to bottom). Blue lines: far-field approximation. 
Red line: exact expression (Fourier-Bessel basis). 
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Fig. 10. Transmittance versus distance normalized w.r.t. far-field limit. m’ = 
2, m = 1 and n = 0...4 (top to bottom). Blue lines: far field approximation. 
Red line: exact expression (Zernike basis). 


impact the decay rate but only modifies an amplitude factor. 
For n < m, the amplitude is proportional to m~° in the 
far-field. So the question arises about which decomposition 
should be chosen for the fields over the apertures. Regarding 
the far-field behavior we may prefer the Zernike approach, 
which allows one to rapidly neglect contributions related to 
fast radial variation of fields, as the distance is increased. 


APPENDIX A 


The Bessel functions orthogonality is defined on the unit 
disk as: 


Onl dys Aa ai (Ar) 


(21) 


1 
l Fm,n(Po) Fi npo) po dpo = 
O 
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with a closed form spectrum given by: 


~ 


1 
F mab) =| Fm,n(po)Jn(6po) podpo = 


22 
— XAnmIn—=1(Anm) In (8) 
Apna)” = p 
APPENDIX B 


Starting from the asymptotic approximation of the Bessel 
functions (valid for high value of the argument): 


(23) 


The m-th positive zero of J,,(x) can be approximated (for high 
values of m) as: 

T T 
Anm S zn +1) +7 +m- i)i 


After inserting (24) in the asymptotic form of the Bessel 
function (23), we get: 


(24) 


Jail Anm) ~ 





COs (S(n+1)+ Í + (m -— 1)r 








T Anim 2 
(n-l)t ~r 
a F) 
2 
A cos (7 + (m — 1)r) 
2 mM 
E TAwm io) 


(25) 


Now, the dependence of the transmittance on m (expression 
between brackets in (20)) can be written as: 


— 2 1 Layee" 


no T (Anra ami e 


Finally, using the approximation (24) of Anm, for high values 
of m in comparison with n, the dependence on those param- 
eters are m—?/? and (m’/)~3/?, 


(26) 


APPENDIX C 


Let us consider two antennas. The first one as the transmitter 
and the second one, the receiver (see Fig. 11). Referring to the 





(a) (b) 


Fig. 11. Monoporte modeling of two antennas. a): the transmitting antenna 
with impedance Zr connected to a voltage source V, with Zg being the 
source impedance. b): the receiving antenna with impedance ZR on which 
an impinging electric-field modeled with a voltage source Voc, and Zz, being 
the load impedance. 


Poynting theorem, the transmitted power is defined as: 


1 x 
Pr= zlte {VrIp} (27) 
where Re stands for the real part, Ir = V/(Za + Zr) and Vr 
are respectively the current and the voltage at the transmitted 
antenna port. Expression (27) can then be rewritten as: 


1 vi 


Pr = ~R. {Zr} [Zr + Zale 


— 28 
5 (28) 
The available power P,,, 1.e. the power delivered by the 
source can be obtained from (28), after inserting the matching 
condition Zp = ZG, which leads to: 


1 VE 


To~ = 


1 |V? _ A Mrllér + Ze 
8 Re {Za} 8 


Re {ZG} 


The current |/7| at the transmitted antenna port can then be 
written as 


(29) 


Re {Za} 


In|? = 8P 5 
Irl Zr + Ze|? 


(30) 


Similarly, the received power on the load is: 


l [Voc] 1 Re {ZL} res 2 
Pee ee ee i 
=~? (413 Zrt 2r t2 | rl? m 
Expression |V,.|?/|I7|? between the brackets is the square 
of the absolute value of the transmittance |H |? defined in this 
paper. After replacing expression of |Ir|? given in (30), we 
can then rewrite (31) as 


Re{Zry. Re{Ze} 
Pa =4 4 P, —— a 32 
R | | Zr + Z|? Zr + Za)? ( ) 
Thereafter, the following identities can be easily proven: 

Re {ZL} 7 1— |['R|? (33) 

Za + ZL? AR. {ZR} 

R.{Z 1—|Cr/? 
{Ze} _ 1- [Pr ae 


[Zr + Ze 4B. {Zr} 


where Dp = (Z, — Zh) /(Zr+ Zz) and Pr = (Zr — 
Ze:)/(Zr + Zq@) are the reflection coefficients. 

Inserting (33) and (34) into (32), the received power Pp 
can be related to the available power Pav as follows: 


1 


- 1 ypPp,.a— |r) — el’) — 
Pr = HP Pa = PrP) - Wal) pape {Za} 


This expression allows us to link the received and the available 
power with the transmittance given in (1). 
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